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We study a generic reaction-diffusion model for single-species population dynamics that includes 
reproduction, death, and competition. The population is assumed to be confined in a refuge beyond 
which conditions are so harsh that they lead to certain extinction. Standard continuum mean field 
models in one dimension yield a critical refuge length Lc such that a population in a refuge larger 
than this is assured survival. Herein we extend the model to take into account the discreteness 
and finiteness of the population, which leads us to a stochastic description. We present a particular 
critical criterion for likely extinction, namely, that the standard deviation of the population be equal 
to the mean. According to this criterion, we find that while survival can no longer be guaranteed for 
any refuge size, for sufficiently weak competition one can make the refuge large enough (certainly 
larger than Lc) to cause extinction to be unlikely. However, beyond a certain value of the competition 
rate parameter it is no longer possible to escape a likelihood of extinction even in an infinite refuge. 
These unavoidable fluctuations therefore have a severe impact on refuge design issues. 

PACS numbers: 05.65.-|-b, 87.23.Cc, 05.40.-a 

I. INTRODUCTION 

Population dynamics is a venerable and important subject that has been studied thoroughly for decades Un- 
derstanding ecological systems is interesting in itself, as are applications of this understanding to a large variety of 
important practical problems such as, for example, the spread of a virus or other disease over a host population 01, 
and strategies for the elimination of pests or for the protection of endangered species These examples point to 
the importance of the study of the extinction of populations and of the conditions that lead to extinction as well as 
those that protect against it. 

The oldest population models pose relatively simple rate equations for a total population, e.g., those of the logistic 
variety, and typically take basic events such as births and deaths into account through appropriate rate processes. A 
generic example is the logistic form 15| 

where X is a measure of the size of a population and is therefore nonnegative. The parameter n is the difference 
between birth and death rates of the population, and X = 0, fj, are its steady states. All that is required for survival 
is that ^ > 0. Fairly early on, the possible role of fluctuations in these models was recognized and included by 
way of additional fluctuating terms. The sources of the fluctuations included in this way in the earliest models were 
associated with variations in the external environment in which the population evolves, and appeared as a fluctuation 
in the parameters of the equations. For example, a fluctuating birth and/or death rate in the logistic form might lead 
to the description Q 

^ = ^,X-X' + mX, (2) 

where £,{t) is a zero-centered Gaussian white noise. The role of "internal fluctuations" arising from the fact that the 
populations are finite and discrete was also recognized, and the stochastic contributions to the population equations 
were obtained from an appropriate expansion of a master equation i. Writing the logistic equation in appropriately 
rescaled form, one would have with Y = X/fi and r = /ii, 

dY 



Yil-Y)+a^Y{l-Y)at). (3) 

The parameter a is also determined by the rates that determine fi. Other sources of stochasticity have also been 
recognized, namely, those arising from natural catastrophes and from changes in gene frequencies |a| ■ 
The per capita growth rate of the population according to the deterministic logistic model is 



(4) 
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which clearly decreases to zero as X approaches the steady state value /i but that has a maximum as X ^ 0. A point 
that was observed early on in the biological ecology literature is that this latter behavior is not descriptive of real 
populations in many cases, and that in fact the per capita growth rate also tends to decrease when the population 
becomes very small 0,0 la- Various possible mechanisms have been discussed for this in the literature, but many of 
them are based on the recognition that not only the total population but also its density may play a crucial role in 
population dynamics, particularly in questions related to extinction. 

More recent population models have recognized the importance of describing not only the total populations but also 
its spatial distribution, thereby getting a handle on density-related issues. Perhaps the easiest way of dealing with 
a spatial distribution is by wa y o f reaction-diffusion equations, and indeed such studies in the biological/ecological 
context have a long history [3, IHL IH, Hj, HI, Hll • Reaction-diffusion equations describe deterministic continuous 
densities. Therefore, whereas they do deal with some of the density-dependence issues of great importance in the 
assessment of extinction probabilities, they do not in general include the other feature discussed above that is also 
of great consequence in these problems, namely, the effect of fluctuations. These fluctuations can again arise from a 
variety of internal and external sources. Herein we address the problem of internal fluctuations. 

Populations are composed of discrete numbers of individuals. That populations are not continuous is particularly 
relevant at low densities and at low total populations. The importance of these effects have been recognized and 
analyzed in various ways in recent years, most prominently in works that deal with a master equation description of 
the population instead of a reaction-diffusion model |la . Il7| . These master equations are typically solved numerically 
and provide important information on various measures for extinction. 

We take a different approach to the problem, our interest lying in an analytic formulation of the problem of internal 
fluctuations and their consequences. On the one hand, an analytic approach necessarily requires a more narrow focus 
on a particular question, and therefore requires some restriction as to the range of issues that we can address. On the 
other hand, analytic results may provide insights that are difficult to obtain in purely numerical work. 

Our specific question arises from the following situation. Suppose that a population lives in a one-dimensional 
bounded refuge with favorable conditions for life, while outside of the refuge the conditions are extremely harsh. A 
question that has been considered prominently in the mathematical literature is that of the critical size of such 
a refuge that will guarantee the survival of the population [Tol IT5|. that is, the critical size beyond which the 
extinction probability vanishes, Pext = 0. This problem, and many generalizations of it, have been solved, but to 
our knowledge existing analytic solutions are based on the assumption of a continuous population. We show that 
accounting for the discreteness and finiteness of the population profoundly modifies the known results and may reopen 
an array of questions for reconsideration. 

Our starting point is a master equation to describe the following scenario. In a refuge of volume V we have a 
population composed of a single stream of individuals called A that diffuse with constant diffusion coefficient D on 
a d-dimensional lattice (later d is set to unity) with lattice spacing Ax. We further suppose that the individuals 
reproduce by fission {A ^ A + A, as may be the case with bacteria) with rate constant cti, they die {A 0) with 
rate coefficient a2, and they compete according to the reaction A + A ^ A with rate coefficient A. We suppose 
that (Ti > (72, because otherwise the population eventually becomes extinct with certainty. We also assume that the 
nutrients are homogeneously distributed so that the coefficients do not depend on the lattice site. The microscopic 
description of the system within the refuge is given by the master equation 

dP{{n,};t) v^f D ^ ^^.^^ ^ 

j {m} 

- niP{. . .,ni,nm, ■■■;t)] 
+ A [(rij -I- l)njP(. ..,n^ + l,...;t)-n^{ni~ 1)P(. . . , n^, . . . ; i)] 
+ CTi [(n^ - 1)P(. . . , - 1, . . . ; i) - niP{. . . , n,, . . . ; i)] 

+ (72 [(n. + 1)P(. . . , + 1, t) - niP{. . . , n„ . . . ; t)] i. (5) 



Here the occupation number nj is the number of bacteria at site j and {rij} =(..., nj, Jij+i, . . .) is the set of all 
occupation numbers. The index i in the sum ranges over all lattice sites and m denotes the set of nearest neighbors of 
i. We choose a homogeneous initial condition given by an uncorrelated Poisson distribution with average population 
n at each site: 

P(K};0) = e-"n— ■ (6) 

I 

Outside of the refuge living conditions are assumed to be extremely harsh, represented by the addition of the process 
A ^ with a large rate coefficient 7. In the limit 7 ^ ex: (certain death outside of the refuge), we can implement this 
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contribution by requiring that P{{ni};t) vanish whenever Uk 0, where k denotes any site just beyond the refuge 
boundary. Thus, "harshness" is implemented as a set of boundary conditions. 

In Sec. mwe discuss the derivation of a stochastic Langevin equation whose moments can be related to those of 
the fluctuating population described by the master equation. In Sec. lIIII we present criteria for survival or extinction. 
In this section we discuss whether a refuge can be made sufhciently large to make extinction unlikely. A summary of 
our results and some directions for future research are presented in Sec. IIVI 



II. THE MODEL 



It is well known that the continuum (mean field) description of the local concentration p{r,t) is the Fisher- 
Kolmogorov-Petrovsky-Piscunov (FKPP) equation [l8| 



dpir,t) 
dt 



DV^p{r,t) + {ai - a2)p{r,t) - \p{r,tf, 



(7) 



which has been the phenomenological starting point in a huge number of problems in physics, chemistry, and biology. 
In one dimension and with the refuge extending inside the interval the equation is to be solved subject to 

the boundary conditions 



L 



0. 



(8) 



For small L the only solution to this problem as t — > cx) is p{x, t) — > 0, that is, the population becomes extinct. There 
is a critical refuge size Lc beyond which extinction does not occur and the population will be nonzero (albeit small). 
Since p{x,t) is small near to find this critical value it is appropriate to linearize Eq. ((TJ around zero population 
density and solve the simpler equation [T^ 



dp{x,t) _ ^d^p{x,t) 



+ (cTi - (T2)/o(a;,t), 



dt dx^ 

subject to the boundary conditions The Fourier decomposition of the solution is given by 

p{x,t) = y^an/9„(a;,t). 



(9) 



p„{x,t) = exp 



-D t 



sm 



riTT 

T 



(10) 



where a = ai — (72, n = 1, 2, and the a„ depend on the initial condition. Thus, all Fourier modes decrease to zero 
in time if and only if L < 7ry^, leading to the conclusion that for these values of L, p„(x,t) when t ^ oo for 
any n. We can therefore conclude that the critical length of the refuge is 




(11) 



That is, the population becomes extinct with certainty ii L < Lc and it certainly does not become extinct ii L > Lc 
(in fact, the Fourier components grow without bound in this case, but the linearized equation is then no longer valid). 
This result, and extensions of it, have been known for more than fifty years j^, and has served as a background for the 
design and analysis of refuges. Note that the linearization argument that leads to Eq. ^ is valid for any nonlinearity 
that can be neglected near extinction in Eq. (0, and that therefore the resulting Lc is obtained not only for this 
specific model il4i] . 

The above mean field analysis does not take into account the effect of the fiuctuations associated with the fact that 
the population is really discrete and finite. To do this short of solving the full master equation, one must formulate 
a generalization of the FKPP model that takes the resulting fluctuations into account. Such gen eralizations exist in 
the literature in other contexts. In particular, on the basis of a theory first proposed by Doi [l^, further elucidated 
by Peliti [2^, used for the study of critical phenomena associated with bulk transitions in reaction dynamics [2ll |. 
explicitly applied to a reaction-diffusion front problem by Pechekin and Levine p^ . and subsequently used in a 
number of other contexts related to Fisher waves (see e.g. H^), one arrives at a stochastic differential equation for 
a field whose moments can be related to those of the population in the original master equation. The connection 
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comes about as follows. The master equation Eq. lO can be projected onto a problem in quantum field theory by the 
construction of an adequate Fock space. Since this calculation has been performed many times in the literature we 
only indicate some of the main steps. A detailed review of the procedure can be found in [23|. Consider the following 
operator algebra: 



,a+] = (5„-, [ai,aj]^0, [af,aj]=0, 



(12) 



where at, af are destruction and creation operators and the square brackets denote the commutator [20j . A state 
vector can be defined as 



|V'(t)) = ^P(ni,n2,...;t)ar^4"^-|0) 



(13) 



where P(ni,n2, ...;i) = P({n,;},i) is the solution of the master equation, and the sum is performed over all possible 
configurations of the {rii}. The state vector obeys the imaginary time Schrodinger equation 



whose formal solution is 

In our case the (non-hcrmitian) Hamiltonian is 

D 



d\m) 

dt 



= -Hm)) 



\m) 



(Aa;)2 



aliam - ai) - A(l - af)afa^ + ai[l - af]afai - 0-2(1 - al)ai 



{m} 



Note that we recover Eq. Q by substituting Eqs. ifT^ and into Eq. (fTHl . 
One introduces the Glauber state as the projection state 



N 



(14) 



(15) 



(16) 



(17) 



The expected values of observables A{{ni}) can then be written as 

(A(t))^^A({n,;})P({M;t) = (5|A 



-Ht 



\m) 



where A is the operator obtained by replacing every rii in the function A{{ni}) by the number operator fii 



(18) 



The steps that we do not repeat |21|, |22| show that this second-quantized theory can be expressed in terms of a 
path integral. Furthermore, from this path integral one can derive a "classical" action from which one can in turn 
derive the mean field continuum equation Eq. Q. If one concentrates on the path integral itself instead of moving 
on to the classical action, one can show an equivalence between the master equation and the following Langevin 
equation pll2^: 



dt 



+ v/2 [aiij{r,t)-X^^r,tmr, t), 
where ^{r,t) is Gaussian white noise with mean and correlation given by |24j : 



(19) 
(20) 



(e(r,i)C(r',t')) =<5(i-mr-r'). 
The multiplicative noise in Eq. (|19|l must be interpreted according to Ito. 



(21) 
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It is important to note that the noise in Eq. H19(l can be imaginary, so the value of ip is in general complex. 
This makes it clear that the field is not to be interpreted as the population density. The field ip is the complex 
eigenvalue of the destruction operator, so one can relate its moments to the moments of the population density using 
the commutation relations Eq. (|12|l and the property of the Glauber state (5*1 a+ ~ {S\. For the first two moments 
we have: 

(p(r, t)) = (n,> = (a+a,) = {S\ a+a, = (S\ a, = (a,) = (V(r, t)) , 

(p^(r, t)) = {n^ni) = (alaiala^) = {alajaiai + a+Oi) = {S\ af^af + a^ai \^p{t)) 

= {S\ a\ m)) + {S\ a, m)) = (a^) + (a,) = {i,\r, t)) + (V'(r, t)) . (22) 

Higher moments can also be calculated following this procedure. Here the lattice site i is associated with the volume 
(Aa;)"^ around point r and we have set Ace = 1 for economy of notation. 
As a side note we point out that the so-called stochastic FKPP equation of the form 



dp{x, t) 



can not represent the population density either since it does not allow fluctuations around the equilibrium population 
density p = ci/A because the population is real, while these fluctuations do not vanish in the original discrete 
system. However, a connection between Eq. H23|l and the original master equation Eq. @ without the death process 
has recently been elaborated 0|. This connection is different from the one presented above involving the stochastic 
equation l(Tn|l. 



III. EXTINCTION PROBABILITY 



The mean fleld model for the refuge population predicts certain extinction if the refuge size L < Lc and certain 
survival if L > Lc- In this section we show that the fluctuations destroy this certainty of survival, that is, that survival 
is in fact never guaranteed. 

In the biological literature one encounters the concept of a minimum viable population |^ IT^ . There are some 
variations in its definition, but one that seems to be widely accepted is that a minimum viable population is the 
smallest isolated population having a 99% chance of remaining extant for 1000 years despite the presence of foreseeable 
fluctuations. Clearly one can vary the two numbers that appear in this criterion, but, in any case, it involves a 
probability of survival in a finite period of time. We wish to make clear at the outset that our analytic theory is not 
yet developed to the point of being able to handle finite times, and we are therefore not able to implement a criterion 
of this sort. Instead, the theory (like the mean field prediction stated above) deals only with steady state probabilities. 

Near extinction one can show directly from the master equation @ that the competition term (the term proportional 
to A) can be neglected. As a result, in this limit the contributions A^^ in the stochastic equation (|19|l can also be 
neglected, and to find the critical refuge size one only needs to deal with the simpler Langevin equation 

^^^^ = DV^^ix, t) + a^{x, t) + v/2aiV(a;,i)e(x, t). (24) 

The mean value of tp obeys the equation 

dmx,t)) ^ dm,{x,t)) , 

Ql = ^ + 'j{Hx,t)), (25) 

which is the same equation, Eq. Q, that we solved for the deterministic case. Since {ip(x, t)) = {p{x, t)), we have the 
same boundary conditions as in Eq.©, which implies that the mean value of the population density is zero, that is, 
that extinction is certain, for L < Lc- Since {ip(x,t)) — {p{x,t)) > for L > Lc, this in turn implies that extinction 
is certain if and only if L < Lc- This, however, does not tell us whether or not the extinction probability is zero for 
L > Lc (as predicted in mean field theory). Indeed, we now show that the extinction probability is greater than zero 
for any finite length of the refuge. 

It is fairly obvious that Pext does not vanish if one has a finite population and thus a nonzero probability (albeit 
perhaps small) of a total extinction event. We will place this statement on a more formal footing. However, an 
arbitrarily small extinction probability is not necessarily very important since in most design problems one deals 
with at least some uncertainties. It is therefore useful to establish a criterion as to what constitutes a non-negligible 
extinction probability and what the minimal size L* of the refuge must be to insure a probability smaller than this. 
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Clearly, L* > Lc- We will establish such a criterion and show that (depending on parameter values) the size L* may 
be much larger than and even infinite. 

First, to establish the obvious fact that fluctuations lead to a nonzero extinction probability even for L > Lc, let 
us suppose that we begin with an infinite refuge and exclude death events ^ — > altogether. The population inside 
an interval of length L of this infinite refuge is then greater than or equal to the population inside our actual finite 
refuge of length L. The population density inside the finite interval of the infinite refuge in the steady state is given 
by the Poisson distribution 

P({nJ) = e--An^^^^' (26) 

SO that the probability of an unbounded realization of the process is zero. Since the integral of a bounded function 
over a bounded interval is bounded, we can conclude that the total population inside the interval is finite. Therefore 
so is the total population inside the finite refuge. Furthermore, with ^-correlated Gaussian fluctuations ll^H^, even 
the most rare fluctuations occur in a finite mean time and therefore a rare fluctuation caused by the reaction A fl) 
could kill the entire population in a finite time. We thus conclude that the probability of extinction is greater than 
zero whenever L < oo. In the case of an unbounded refuge the mean value of the population density is cr/A and is 
homogeneously distributed, so a local fluctuation cannot kill the population. We can summarize these conclusions as 
follows: 

Pext = 1 if L < Lc 

< Pext < 1 if Lc < L <oo 

Pext when L -> oo. (27) 

We can furthermore state that since the number of individuals in the population increases as the length of the refuge 
increases, the probability of extinction decreases monotonically and continuously as L increases continuously. Thus, 
consideration of the internal fluctuations drastically changes the picture; now there is no finite size of the refuge that 
can be considered absolutely safe for the population. However, we have not been able to derive the explicit functional 
dependence of Pext on the size of the refuge. 

As an alternative to such an explicit full solution, we propose a criterion for deciding when a population is under 
substantial risk of extinction, and calculate the critical size L* of the refuge associated with this criterion. Suppose 
that a population density p{x) in a refuge of length L described according to mean field Eq. Q has a maximum 
steady state density AI > 0. Clearly, L > Lc- The maximum occurs at a; = 0, in the middle of the refuge, that is, 
p{x — 0) — M Now consider the description of this population that includes the fluctuations, and compare 

the average population density {tpix)) — {p{x)) with the deterministic result p{x). This average distribution also has 
a maximum at a; = (see below), {il}{x — 0)) = M' . We define the critical size L* as the refuge length associated 
with two conditions: 

1. The deterministic and stochastically obtained maxima are equal, Af — M'; 

2. The standard deviation of the stochastic population density is equal to its mean at each x. 

Clearly, our idea is that a population density whose standard deviation everywhere equals its mean is in danger of 
extinction, and we call this a "critical population." It is not clear a priori that these conditions can be simultaneously 
satisfied, but in fact it turns out that they can. The second condition above says that for the critical population 

{p{x,t)f ^ {p\x,t)) ~ {p{x,t))\ (28) 

or, equivalently, in terms of the field "0, one for which 

{i;{x, t)f = {i:\x, t)) + (V^x, t)) - (^(x, t)f . (29) 

From Eq. H19|l we obtain the following equation linking the first and second moments of Tp: 

^i^tlH ^ D^^{i,{x,t))+a{^{x,t)) - \{i^{x,tf). (30) 

Inserting Eq. (|29|l for the critical field then results in a closed equation for the first moment, 

= W2(V(x,i)) + (a + m{x,t)) - 2A(0(x,t))2. (31) 
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Note that the apparent irrelevance of the fluctuations in Eq. (|19|l other than their mean value in arriving at Eq. H31() 
is illusory since the relation H22I) and consequently Eq. (|29|l are intimately connected to the precise form of the 
fluctuations. 

In the steady state Eq. H31|l reduces to the boundary value problem 



D 



+ {n + X){i,{x))^2\{4,{x))^ ^Q, 



(32) 



with boundary conditions as in Eq.Q. To find L* we multiply Eq. (|32() by d{'ip{x)) /dx and rewrite the result as 

2 



from which it follows that 



d 
dx 



D 



dx 



{4,{x))' 



2A 



0, 



( d{^{x)) 
2 \ dx 



A) 



{i,{x)y^ 



2A 



constant. 



(33) 



(34) 



Symmetry considerations show that the solution of this problem has only one maximum at a; = '27] , and since the 
first derivative of the solution should vanish at any maximum, we can explicitly write the constant as 



D 



2 V dx 



a + X 



{i^{x)f 



2X 



mx)r 



a + X 



2X 



(35) 



where M is the maximum introduced earlier. Since this maximum occurs at a; = 0, we can integrate Eq. (|35|l for 
x > to obtain 



M 



dz 



(36) 



Taking into account that the population vanishes at the edge of the refuge we can obtain the value of the length L of 
the refuge as a function of the maximum of the average population density: 



L = ^2D 



M 







(37) 



With the change of variables w — z/M and setting e = X/cr we can rewrite this relation as 

r-l 



-L : 



1 



1 



- e 



dw 



1 



(38) 



The steady state deterministic equation obtained by setting the left hand side of Eq. equal to zero is formally 
identical to the critical field equation Eq. H32|) with modified parameters. Since we are choosing L* as the value where 
the maximum of the critical average population density is equal to the maximum of the mean field population density, 
the mean field model provides us with a second relation between M and L: 



D V e 



dw 



y^l - w2 - |e(l - w3)M 



(39) 



The simultaneous solution of Eqs. (|38|l and H39() then leads to the value of L* (and incidentally also of M). 

The two equations reduce to the same equation if A = (e = 0), and in this case L* = Lc and M diverges . The 
reason for the equality is that when the nonlinearity is not present, the total steady state population is unbounded 
(thermodynamic limit) and there are no fluctuations. In general, the two simultaneous equations can only be solved 
numerically, and the results are shown in Fig^ 

A number of important points should be noted about our results. First, L* is indeed greater than Lc except for the 
single trivially coincident point when A = 0, that is, when there is no competition. Thus, according to our criterion 
whereby the risk of extinction is substantial when the mean and standard deviation of the population density are 
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FIG. 1: Left scale: \J a jDh^ (dotted line) and ^ a jDL* (solid line) vs A/cr. Note that ^Jcr/DLc = tt for all values of A/cr, 
while L* grows monotonically up to A/cr = 1, where it diverges logarithmically. Note also that Lc — L* when A = 0, and that 
this is the only common point of the lengths. Right scale and dashed curve: maximum M of the average population density as 
defined in the text vs A/cr. 

equal, one must increase the size of the refuge considerably beyond the deterministic refuge size to avoid this risk. 
Indeed, when A/cr = 1 the size L* diverges logarithmically and the population density becomes uniform at density 
ct/A. As a grows beyond a it is no longer possible to avoid considerable risk of extinction for any finite size of the 
refuge. The increase in Af with decreasing L* can be understood by moving from right to left in the figure: decreasing 
critical refuge size is associated with the requirement of a larger population (and hence a larger maximum population 
density) to insure against extinction. At Lc the maximum M diverges because no population, no matter how large, 
can be protected against extinction below this length. We can obtain an analytic expression for L* for small values 
of A/cr: 



We can also see analytically that L* must diverge when A/ct = 1: at this point the two integrals for L can only be 
equal if they both vanish (which they do not) or if they both diverge (which they do when Af — 1). 



The main point of this paper has been to assess the effects of the inevitable internal fiuctuations on the prediction 
of the risk of extinction of a population in a refuge as a function of the size of the refuge. In the usual mean field 
deterministic model in one dimension, the population is treated as a continuum and one obtains a critical refuge 
length Lc such that the population becomes extinct with certainty if the length of the refuge is below this critical 
value, L < Lc, whereas survival is certain if L > Lc- It is well known that, as predicted by the mean field approach, 
the smaller the area occupied by a population, the lower are its chances for survival |23. However, it is also known 
that a small size or low density of a population may also increase its chances of extinction, a feature that mean field 
theories do not capture 3J . Our approach deals with both of these features simultaneously through the interplay of 
discreteness and nonlinearity in a spatially distributed model. We have argued that the discreteness of the population 
and its finite size make fluctuations unavoidable, and have shown that whereas L < Lc still guarantees extinction, 
L > Lchy no means guarantees survival. While we have not been able to find an analytic expression for the survival 
probability as a function of L, we have focused on a sensible risk criterion for which we have found explicit results up 
to quadrature. Comparing deterministic and average stochastic population densities with a given maximum M, we 




(40) 



IV. CONCLUSIONS 
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have chosen to define a new critical length L* as one for which the mean of the population density and its standard 
deviation are equal, reasoning that this variability implies considerable susceptibility of the population to extinction. 
We have shown that L* > and, most importantly, that L* diverges when the competition rate coefficient grows 
beyond that of the net growth rate of the population. For some parameter values it is possible to protect (albeit not 
with certainty) a population from extinction by placing it in a sufficiently large refuge, one that is certainly larger 
than that predicted by the standard mean field model. For other parameter values (in particular, when competition 
is too strong), it is not possible to evade the risk of extinction (at least according to our criterion) by enlarging the 
size of the refuge. Clearly, these results have serious implications for the expectations in the design of refuges. 

Many questions still remain to be answered. For instance, an exact result for the extinction probability would 
clarify many issues. So would the ability to obtain time dependent solutions so as to deal with a more realistic 
criterion of survival over a long but finite time interval. Also, we have only considered the simplest most generic 
situation, whereas real systems (particiilarly ones not designed in the laboratory) arc likely to be seriously affected by 
many complicating factors, some connected with the necessary extension to higher dimensions. Examples include the 
presence of convection, different boundary conditions, and spatial inhomogeneities inside the refuge. Similar questions 
can and should be posed when multiple species are present. In any case, in light of our results it would seem prudent 
to reconsider other critical size population problems to assess the effects of discreteness. As we have shown in the 
simple model considered here, the consequences can indeed be profound. 
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